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Abstract — One of the main challenges in the emerging smart 
grid is the integration of renewable energy resources (RER). 
The latter introduces both intermittency and uncertainty into 
the grid, both of which can affect the underlying energy market. 
An interesting concept that is being explored for mitigating the 
integration cost of RERs is Demand Response. Implemented as 
a time-varying electricity price in real-time, Demand Response 
has a direct impact on the underlying energy market as well. 
Beginning with an overall model of the major market participants 
together with the constraints of transmission and generation, we 
analyze the energy market in this paper and derive conditions 
for global maximum using standard KKT criteria. The effect 
of uncertainties in the RER on the market equilibrium is then 
quantified, with and without real-time pricing. Perturbation 
analysis methods are used to compare the equilibria in the 
nominal and perturbed markets. These markets are also analyzed 
using a game-theoretic point of view. Sufficient conditions are 
derived for the existence of a unique Pure Nash Equilibrium in 
the nominal market. The perturbed market is analyzed using the 
concept of closeness of two strategic games and the equilibria 
of close games. This analysis is used to quantify the effect of 
uncertainty of RERs and its possible mitigation using Demand 
Response. Finally numerical studies are reported using an IEEE 
30-bus to validate the theoretical results. 



Nomenclature 

8 n Set of indices of generating units at node n 
d n Set of indices of demands at node n 
Q, n Set of indices of nodes connected to node n 
D q Set of indices of Consumers {1,2,..., Np} 
Gf Set of indices of generating units {1,2,..., No} 
G w Set of indices of renewable generating units 

{1,2,..., N w } 
N Number of all buses 
N D Number of demands 
N t Number of transmission lines 
Nq Number of generating units 
Npj Number of blocks demanded by demand j 
Noi Number of blocks bid by generating unit i 
N Gl Number of blocks bid by renewable generating unit I 
7„ m Dual variable associated with the transmission capac- 
ity constraint of line n — m 
ipjk Dual variable associated with the maximum capacity 
limit for block k of demand j 
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p n Locational marginal price corresponding to the gen- 
erating unit i or the demand j that is located at node 

n 

a>i Dual variable associated with the maximum capacity 

constraint of generating unit i 
(f>ib Dual variable associated with the maximum capacity 

limit for block b of generating unit i 
<jj Dual variable associated with the minimum demand 

constraint of demand j 
A^,. fe Marginal utility associated with block k of demand j 
^Djk Price bid by demand j to buy power block k 
Xq. Price bid by generating unit i to sell power block b 
Xq. Marginal operating cost associated with block b of 

generator i 

^Djk Maximum power demanded in block k of demand j 
Pjjjk Minimum power supplied to the demand j 
PDjk Power block k that the demand j is consuming 
P&b Power block b that the generating unit i is producing 
Pgu> x Maximum power output in block b of generating unit 
i 

pmax Maximum power output of generating unit i 
p™« Minimum power output of generating unit i 
Pq W Power block b that the renewable generating unit I is 
producing 

Pqw Maximum wind power estimation for block b of 
renewable generating unit I 

Pqh, Minimum wind power estimation for block b of 
renewable generating unit I 

Pcib ■ Mean wind power for block b of renewable generat- 
ing unit I 

Pmn X Transmission capacity limit of line n — m 

Ap„(j) Change of LMP for demand j 

APojk Change of consumption by demand j 

A Gib Uncertainty in block b of wind generating unit i 

KDjk Curtailment Factor of block k by demand j 

S n voltage angle of bus n 

B nm susceptance of line n — m 

I. Introduction 

The goal in a Smart Grid is to integrate large-scale renew- 
able generation and emerging storage technologies, control 
signals to loads to match supply, and dramatically improve 
energy efficiency. Electricity markets, the entity that carries out 
power balance by scheduling power using bids from various 
generating companies, are crucial components that can con- 
tribute to an efficient smart grid. Electricity market models that 
accurately represent the market dynamics are exceedingly im- 
portant. These models must capture the behavior of the dom- 
inant market players such as generators, consumers, and ISO, 



market mechanisms such as gaming behaviors (ex. bidding 
strategies), real-time prices via Demand Response, diverse 
dynamic drivers (e.g., weather, load, fuel prices, and wind 
supply), and physical constraints (e.g., ramping, transmission 
congestion) (TJ. These market models and analysis of the 
market equilibrium directly help to identify various sources of 
price volatility and dependence and sensitivity to uncertainties 
and intermittencies in renewable energy sources and variations 
in real-time prices due to Demand Response. Our focus in 
this paper is on an electricity market model that captures the 
dynamics of market players, physical constraints, perturbations 
in renewable energy resources (RER), and variations due to 
demand response. 

Various methods have been proposed in the literature to 
determine market models ifTO) . ifTTI . flBl . Given the primary 
purpose of market balance, these methods are primarily fo- 
cused on the analysis of market equilibrium. Methods in |10|, 
IfTTI model the electricity market participants subject to spot 
market equilibrium and their own constraints, leading to a 
method termed Equilibrium Programming with Equilibrium 
Constraints (EPECs). In such models, each agent solves a 
Mathematical Program with Equilibrium Constraints (MPEC) 
02), and by using stationarity theory for MPECs a standard 
mixed complementarity problem (CP) is shown to lead to 
equilibrium lfl6l . 

Another approach to market modeling employs a varia- 
tional and complementarity Problem (CP) formulation (see 
IfTTI . Ifl8l . 1521 . Ifl9l ) which models the equilibrium of the 
electricity market where forward and spot decisions are made 
simultaneously. These models may also be specified by the 
ISOs objective, maximization of social welfare ifTOl . IfTTI or 
maximization of wheeling or transmission revenue |17|, |19|. 
A single-settlement framework based on a linear complemen- 
tarity problem (LCP) is proposed in [14|, which leads to a 
strategic game in which players have coupled constraints and 
the equilibrium properties are analyzed. Yet another approach 
used for electricity markets is the use of Cournot models. 
Here, the strategies of the generator companies are allowed 
to depend on other market participants, who also act so as 
to maximize their own profits. This in turn allows a more 
transparent dependence of costs of generators and consumers 
profit on the market equilibrium (see GUI . GT) . 11521 . 1221 . 

ma, en, id, a, na, s). 

In the context of a market for smart grid, the question that 
needs to be addressed is how the effect of RERs can be best 
captured in a market model. The integration of RER poses the 
challenge of intermittency and uncertainty, both of which in 
turn can affect the economic planning and operation of the 
overall power system. The uncertainties can impose costs due 
to the limited-dispatchability of intermittent generation |31|, 
the variability in resource availability, and errors in forecasting 
resource availability and loads |34|, lf35l . l36l . There is a 
considerable literature covering many important aspects of 
wind power ranging from comprehensive integration studies, 
forecasting methods, and technology issues (see O, [32|, 
[33 1, [34), (35)). In the specific context of the incorporation 
of renewable energy in a market, the results of l37l . [38 1, 
I, 140) . BT) . If42l are more pertinent, all of which deal 



with a perturbed market due to RER. In these papers, a 
stochastic framework is used to capture both overestimation 
and underestimation of available wind power on the optimal 
expected profit of wind power producer, the overall market 
efficiency, and the overall operation cost. In this paper, we 
utilize a deterministic version of the the model developed in 
fiTll and analyze the perturbed market and quantify the effect 
of wind uncertainty on the market equilibrium. 

In recent years, there is a significant interest in Demand 
Response, a concept that allows the electricity demand to 
response to fluctuations in the electricity price [37 1, [43 1. One 
way to reduce wind integration costs is to introduce demand 
response in the form of time-variant retail electricity rate, such 
as real-time pricing (RTP). RTP can potentially reduce wind 
integration and forecast error costs, since consumer demand 
could be made to follow the supply of wind generation by 
using a price signal. Under RTP, if available wind generation is 
less than forecast, the high cost of deploying ancillary services 
to cover the generation shortfall can reduce electricity demand 
and the cost of serving the load [30]. Similarly, because wind 
generation has zero marginal cost, electricity demand may 
increase when there is more wind resource available than 
forecast. 

In this paper, we discuss the equilibrium of an electric 
energy market in the presence of renewable energy sources 
and demand response. Starting from the framework proposed 
in lfl4l . we use the approach based on linear complementarity 
problem to delineate the underlying market model. Repre- 
senting the effect of forecast errors in renewable energy as 
an uncertain parameter, we evaluate its effect on the market 
equilibrium. In a similar manner, focusing on RTP as the 
means by which demand response is realized, and representing 
the effect of RTP as a parametric curtailment on the demand, 
the effect of this parameter on the overall market equilib- 
rium is also studied. A comprehensive market framework 
is introduced in order to evaluate these uncertainty effects. 
This market framework consists of three main participants, 
generating companies including renewable energy sources, 
consumer companies that are capable of responding to market 
changes, and independent system operators (ISO) who clear 
the market by maximizing social welfare [14] and [44 1. The 
overall market equilibrium is first evaluated in the absence of 
uncertainties due to either forecast errors in renewable energy 
sources or due to demand response, and sufficient conditions 
for the proposed market equilibrium to be identical to the Nash 
equilibrium is established. The market response in the presence 
of uncertainties, and in particular, the market equilibrium is 
analyzed using perturbation analysis. A parametric charac- 
terization of the equilibrium shift is delineated. In both the 
nominal and perturbed markets, the equilibriums are derived 
using KKT conditions We next define formally the closeness 
of two strategic game and relate the equilibrium of close games 
using the notion of a — approximation and e — equilibrium. 
IEEE 30-bus system is used to evaluate the results of both 
uncertainties in renewable energy and in demand curtailment. 

The remainder of this paper is organized as follows: In 
Section |n] we present the necessary preliminaries, Section III 
describes the model of the three market participants. In Section 



|IV| the overall market equilibrium under nominal conditions 
is formulated. A perturbation analysis is introduced in Section 
|V| to address the effects due to uncertainty, and the resulting 
effects are summarized. In Section [VI] numerical simulation 



results are reported. Section VII includes concluding remarks. 



II. Preliminaries 

In this section, we provide some preliminaries related to the 
wholesale energy market structure. These include the funda- 
mental theorems of convex optimization which are presented 



in Sections II-A Given the close relationship between convex 



optimization and Game theory, we present in Section |II-B| 
related definitions and theorems related to game theory and 
Nash equilibrium. In particular, the link between the convex 
optimization problem and uniqueness of the Pure Strategy 
Nash Equilibrium is presented including sufficient conditions 
for the latter. These in turn are directly used in establishing the 
equilibrium of the wholesale market in Section [III] We start 
with a few basic definitions. 

Definition 1: A set K C K is convex if for any two points 
x,y £ K, 



ax + (1 — a)y £ K, Vie, y € K and a € [0, 1]. 



(1) 



Definition 2: Given a convex set K C K and a function 
f{x) : K — > R; / is said to be a convex function on K if, 

Vx,y e K and a £ (0, 1), 



f(ax + (1 - a)y) < af(x) + (1 - a)f(y), 



(2) 



Furthermore, a function f(x) is concave over a convex set if 
and only if the function —f(x) is a convex function over the 
set. 

Definition 3: Given a scalar- valued function f(x) : 1" — > 
M. we use the notation Wf(x) to denote the gradient vector of 
f(x) at point x, i.e., 

T 



V/(x) 



dx n 



(3) 



Definition 4: Given a scalar- valued function f(x) : 
nLi K " 1 * ^ Bwe use the notation V,/(a;) to denote the 
gradient vector of f(x) with respect to x t at point x, i.e., 

T 



V</(aO 



df(x) 
fir 1 ' ■ 



(4) 



A. Convex Optimization 

Consider a generic optimization problem 
Maximize f(x) 
s.t. g n (x) = 0, Vn=l,...,iV 

(5) 

^iJmAW^Cm, Vm = l, ...L 

71=1 

where /(ie) is called the objective function or cost function, 
R is a matrix of constants and c m are constants. We assume 
that f(x) : M" — >• M is a convex function to be maximized 
over the variable x, the functions g n (x) as equality constraints 
are affine, and the functions h n (x) as inequality constraints 



are cocave. With these assumptions the optimization problem 
^ is termed a convex optimization problem. In addition, the 
constraint set for the optimization problem is convex which 
allows us to use the method of Lagrange multipliers and the 
Karush Kuhn Tucker (KKT) theorem which we state below 

ma, nm. 

Theorem 1: Consider the optimization formulated in d3), 
where f(x) is a convex function, g n {x) are affine functions, 
and h n (x) are cocave functions. Let x* be a feasible point, i.e. 
a point that satisfies all the constraints. Suppose there exists 
constants A„ and fi m > such that 



V/(x*) + ^A n V. 9n (x*)- 



n=l 



^ H m (R m n^h n (x*) - Cm) =0 Vn = 1...JV 



(6) 



rn — 1 



V- m {Rmnh n {x*) - C m ) = ,VTO=1, 



then is a global maximum. If f(x) is strictly concave then 
x* is also the unique global maximum. 
Proof: see 



B. Game Theory and NASH Equilibrium 

An alternative game -theoretic approach can be used to 
describe the market equilibrium ||5D . ll52l . To describe a game, 
there are four things to consider: 1) the players, 2) the rules 
of the game, 3) the outcomes and 4) the payoffs and the 
preferences (utility functions) of the players. A player plays 
a game through actions. An action is a choice or election 
that a player takes, according to his (or her) own strategy. 
Since a game sets a framework of strategic interdependence, 
a participant should be able to have enough information about 
its own and other players' past actions. This is called the 
information set. A strategy is a rule that tells the player which 
action(s) it should take, according to its own information set 
at any particular stage of a game. Finally, a payoff function 
expresses the utility that a player obtains given a strategy 
profile for all players. More formally a strategic form game is 
defined as follows. 

Definition 5: A strategic forms game is a triplet 
(X, (Si)i£.z, {ui)i£i) such that I is a finite set of players , 
i.e. X = {1, /}, Si is the set of available actions for player 
i, Si £ Si is an action for player i, and finally Ui : S — > R is 
the payoff function of player i where S — Yli Si is the set of 
all action profiles. 

In addition, we use the notation S-i = [sj]j^i as a vector of 
actions for all players except i, S-i = Ylj^ti Sj as the set of all 
action profiles for all players except i, and finally (si, S-i) £ S 
denoted a strategy profile of the game. Informally, a set of 
strategies is a Nash equilibrium if no player can not do better 
by unilaterally changing his or her strategy. 

Definition 6: A Pure Strategy Nash Equilibrium of a strate- 
gic game (1, (Si) ie x, {ui)iex) is a strategy profile s* £ S such 
that for all i £ 1 



Ui{si,s%) < u^s^s*^) Vs, £ Si 



(7) 



Definition 7: Let us assume that for player i G I, the 
strategy set Si is given by 

Si = {x, G R m \hi(xi) > 0} (8) 

where hi : R™ — > K is a concave function and the set of 
strategy profiles S = Yll—i Si C ^-T I s a convex set. 

The payoff functions (ui,...,uj) are said to be diagonally 
strictly concave for x £ S if for every x* , x £ S, we have 

(x - x*) T Vu(x*) + (x* - x) T Vu{x) > 0. (9) 

In the current context, the optimization problem in |5]l can 
be viewed as a game played by all player i £ I, and so that 
together they can arrive at a Nash equilibrium that satisfies the 
condition |6} in theorem [TJ 

We now present a theorem that links the equilibrium point 
of the convex optimization problem in <|3j to Nash equilibrium 
using the above definitions. 

Theorem 2: Consider a strategic form game 
(X, (Si)igi, (ttj)i£i) with Si as in Definition [7] Assume that 
the symmetric matrix (U(x) + U T (x)) is negative definite 
for all x £ S, where the matrix U(x) is the Jacobian of 
the gradient vector of u(x), i.e. the jth column of U(x) is 
du(x)/dxj, j = 1, 7. Then the game has a unique pure 
strategy Nash equilibrium identical to the optimal solution 
x* £ Yii=i ^-T °f tne following optimization problem 

max Ui(yi,x*_A 

s.t. hi(xi) > 0. (10) 

Proof: From Theorem 6 in [29] we have that if U (x) + 
U T (x) is negative definite for all x £ S, then the payoff 
functions {u\,...,Ui) are diagonally strictly concave. This in 
turn, using Theorem 2 in 11291 , implies that there is a unique 
pure strategy Nash equilibrium that is identical to the optimal 
solution x* £ Yii=i °f me convex optimization problem 
in {TO). ■ 

III. Modeling of Market Agent Behavior 

In this section, we derive the model of the overall electricity 
market. The main components of this market include (i) Gener- 
ating Companies (GenCo) (ii) Consumer Companies (ConCo) 
(iii) Independent System Operator (ISO) [14|. We focus our 
attention in this paper on the wholesale market, and use 
a deterministic framework to represent all components. The 
underlying problem that is to be addressed is the determination 
of the power production of the generating units and the power 
demanded by the consumers such that power is balanced at 
all network nodes and social welfare is maximized while 
meeting capacity limits. This is a fairly complex optimization 
problem since each of the components (i)-(iii) are subjected 
to constraints that may conflict with each other. For instance, 
consumers are interested in maximizing their benefit, and are 
therefore interested in cheap electricity. Availability of such 
cheap generation is limited due to capacity constraints on 
transmission lines. In this case, the ISO may need to force 
more expensive units to keep things in equilibrium. The solu- 
tion of the resulting optimization problem not only determines 



the optimal dispatch but also the optimal Local Marginal 
Price which is determined as the corresponding shadow price. 
In what follows, we model each of the components (i)-(iii) 
together with their constraints and the optimization goals. 

A. Generating Company 

It is assumed that the generating company consists of Nq 
generating units, and that the production of each generating 
unit i, i = 1, . . . , Nq is divided into b power blocks, where 
b = 1, . . . , Na it and is denoted as Pcib- The associated linear 
operating cost is denoted as A§ . The goal of the company 
is to maximize its overall profit, and is stated as 

N Gi 

maximize ^ ^ (p„ (i) - \% %b )P Gib (11) 

i£G f b=l 

where denotes the LMP of unit i at node n in the 

network. The power production is subject to the following 
constraints: 

N Gi 

Y^l',;,:,- J>rr :n,: :: ■ Cj (12a) 

b=l 

P Glb < PSa ■ &b;V* e Gf,b = 1, . . . ,N Gi (12b) 
^G^O^eG/;^!,...,^ (12c) 

where ctj and 4>n, are the corresponding shadow prices. The 
decision variables of this problem are the amounts of power 
Pcib to be generated by each unit i in each block b. While 
p n (i) are variables in the overall OPF problem, they can be 
viewed as fixed values with respect to the optimization of ( 1 1 1. 
Necessary and sufficient conditions for the optimization of 
( [TT| ) subject to ( |12a| > - \\2c\ are enumerated below, using dual 
variables ai and <fin,, which correspond to the KKT conditions 

ma 

Paib( - P„(i) + Ag sb +a t + (j) ib ) = Vi G G f ; (13a) 
b=l,...,N Gl 

cti( - £ P Gl b + Par) - 0; Vi £ G f (13b) 
6=1 

<t>ib( P Gtt? ~ Pdb) =0;VieG f ;b = l,...,N G i (13c) 

Definition 8: GenCo i is said to behave rationally if the as- 
sociated linear operating cost denoted as X G is nondecreasing 
with respect to PG ih - That is 

d(\C ib )/d(P G J > 0, W£G f , b=l,...,N Gl . (14) 

B. Consumer Modeling 

A consumer company (ConCo) is assumed to consist of 
No units, and the demand of each unit j, j = 1, . . . , Nd, 
is divided into several blocks Nrj-, with a block-index k = 
1, . . . ,Nd-. The associated linear utility function is denoted 
as which represents the value of using electricity for 

the consumer. The goal of a ConCo is to maximize the total 
profit while consuming electricity. This profit, for a unit j in 
block k connected to node n, is determined as the difference 



between the utility A^,. and the corresponding LMP p n (j) 



[14|. Assuming that the corresponding power consumed is 
denoted as Pojk, the maximization problem can be posed as 



Maximize ^ ^ (^D jk ~ Pn(j))Pojk 

j£D g fc=l 

As before, this is subjected to the constraints 



P Djk < PES ■ <r s ;Vj < » ■<■■'■■ i N Dj 



E P ^ ^ P n!k 
fe=i 



ipjk;Vj g D q ;k = l, . 



,iV, 



(15) 

(16a) 
(16b) 
(16c) 



where <jj and ipjk are the corresponding dual variables for 
(16a i and (16bi. While the constraint ( 16c I is almost always 
satisfied in current-day markets, this may not be the case in 
micro-grid topologies. The decision variables of this problem 
are PDjk, the amounts of power to be consumed by each 
demand j in each block fc. While p n {j) are variables in the 
overall OPF problem, they can be viewed as fixed values, as 



before, with respect to the optimization of (15 1. Necessary 



and sufficient conditions for the optimization of (|T3J subject 
to ( |16a| i and ( |16b| > are enumerated below, using dual variables 
dj and ipjk, which correspond to the KKT conditions |[T4l 



P D 3 k{Pn(j) 

VjeD q ;k = l 



An 



- °i + V'jfe) = o 



(17a) 



*i ( E P W ~ P Djk) = 0; Vj G D q ; k = 1, 
fc=i 



<Jjk 



-P, 



Djk 



r>rnax\ 
r Djk ) 



= 0;VieD,;fc = l, 



(17b) 
(17c) 



Definition 9: ConCo j is said to behave rationally if the 
associated linear operating cost denoted as A^, is non- 
increasing with respect to Pn Jfc . That is 

«.JMP D .J<0, VjeD q ,k=l,...,N DJ . (18) 



C. Effect of Uncertainties and Demand Response 

The most dominant impact of the introduction of distributed 
energy resources is uncertainties, which can directly alter the 
power generated. This in turn can affect the overall market 
equilibrium as well as the corresponding LMP. We take a 
first step in this direction by introducing uncertainties in the 
decision variables introduced in section A. We first separate 
the family of Pcib into Pq^, i = l,...,nc, and P($ lb , 
I = 1, . . . , nji, where nc denotes the more conventional 
dispatchable generating units, and tir denotes distributed 
energy resources such as those based on wind and solar energy, 
which are non-dispatchable. At each node, both conventional 
generators and wind generators may be committed for pro- 
duction. We assume that the wind GenCo are competitive and 
that they submit their bids to the market as other conventional 
GenCo, and not modeled as a negative demand. It is also 
assumed that the available wind energy is overestimated which 



implies that if the assumed power is not available, power can 
be purchased from an alternate source or that loads can be 
shed. Underestimation implies that surplus power is either 
sold to adjacent utilities, or consumed through fast redispatch 
and automatic gain control, or reduced through reduction of 
conventional generation. This case is however not discussed 
in this paper. 

Using the above discussion, the objective function defined 
in ( fTTj ) is modified as 



N Gl 

maximize (Pn(i) 

/gG„ b=l 



Ag™ 



C r w JA Wlb ))P% lb (19) 



where Pqh, is the wind power from the I wind generator, 
and assumed to lie in the interval [Pqu, , Pi 
is due to wind uncertainty, given by 



Glb„ 



P Glb^Glb 



(20) 



where Pqm i s the mean wind power, and < Aqw < 1. We 
note that if there is underestimation, then Acib needs to be 
negative. 

Xq-w is the marginal cost function for the I th wind generator 
which is very close to zero. Finally, C£, AA Wlh ) is a cost 
incurred by committing specific generators as reserves [42|, 
due to the wind uncertainty A Wlb , and is modeled as a 
quadratic function 



b Wlb A Wlh 



^Wlb_ A 2 

2 Wlb' 



(21) 



We define another quantity x w which represents the percentage 
of wind penetration as 



Gib 



j£D q l^k=l r Djk 



(22) 



where P^jk ' s me maximum power demanded by consumer 



j in block k defined in ( 16a i. We note that the impact of wind 
power on the market equilibrium is much smaller if x w is 
small, i.e., if wind penetration is low, than if x w is large. 

Yet another component that we introduce into the picture 
is Demand Response ll53ll . |9l . As for GenCo, categorizing 
all ConCo units into dispatchable and non-dispatchable ones, 
we consider a dispatchable load Pojk and assume that it 
is affected by Demand Response and is sensitive to price 
changes. This effect is modeled using a control parameter 
KDjk, and denotes the response of the consumers to a change 
in the Real Time Price (RTP) which occurs due to Demand 
Response. 

PDjk — PDjk(^ ~ K Djk) < K D j k < 1 (23) 

PDjk denotes the consumption incorporated with demand 
responsiveness into RTP. A positive KDjk, denotes a decrease 
in the ConCo consumption while a negative KDjk, denotes 
an increase. In this paper, we restrict our attention to positive 
KDjk since our focus is on cases where there is a shortfall in 
the non-dispatchable GenCo, i.e. Ag;j, > 0. It is assumed that 
KDjk is suitably calibrated to represent the effect of RTP on 



the consumer behavior. Noting that this control parameter is 
proportional to the elasticity factor e defined as f9l 

(24) 



P»0) 



and is a measure of the consumption response to RTP changes, 
we denote this control parameter, KDjk as Curtailment Factor. 

In contrast to k-Gkb, which may be unknown, curtailment 
factor KDjk> is assumed to be controllable. Nevertheless, the 
equilibrium as well as the LMPs are altered due to the presence 
of the perturbation parameters Aga, and KDjk- In Section III, 
a similar optimization procedure as above and equivalent KKT 
conditions to ( |13a[ >-( fl3c] > and ( |17a[ )-( [T7c] > are derived, and the 
dependence of the equilibrium and LMPs on the perturbation 
parameters are discussed in detail. 

It should be noted that the effect of RTP is modeled in this 
paper in the form of The inherent assumption here is 
that the nondispatchable ConCo observes the LMP, which is 
the solution of (17ai-(17ci, and suitably adjusts its demand. 



This adjustment is represented through the curtailment factor 
KDjk and therefore the elasticity factor e. Therefore, it may 
be argued that the RTP is assumed to be the same as LMP in 
this paper, and that the effect of RTP is modeled essentially as 
a Demand Response parameter, denoted by KDjk. Further im- 
provements of such a modeling of RTP and demand response 
are the focus of ongoing research. 

D. Independent System Operator ISO modeling 

In addition to GenCo and ConCo, a competitive electricity 
market also includes a coordinator which is independent of 
these profit-maximization entities ll46l . B31 . [50|. Such an 
independent entity is necessary in order to maintain order in 
the overall market. Denoted as Independent System Operator 
(ISO), this entity is responsible for guaranteeing a non- 
discriminating access for all participants to the market, and 
enforces various guidelines such as maintaining transmission 
capacity limits. While in practice, distinctions are made be- 
tween an ISO and what is termed a Market Operator (MO) 
where the latter is responsible for financial management and 
the latter for overall physical supervision, in this paper we 
combine the two roles into the same rubric of an ISO. The 
responsibilities of the ISO are to enforce transmission capacity 
limits, to maintain independence from the market participants, 
to avoid discrimination against the market participants, and to 
promote the efficiency of the market. 

The electricity market that we consider in this paper is 
wholesale and is assumed to function as follows: First, each 
generating company submits the bidding stacks of each of 
its units to the pool. Similarly, each consumer submits the 
bidding stacks of each of its demands to the pool. Then, the 
ISO clears the market using an appropriate market-clearing 
procedure resulting in prices and production and consumption 
schedules. The market-clearing procedure may introduce net- 
work constraints, which model losses [46 1 and line capacity 
limits |49l . In this paper we include network constraints in 
the underlying market model and the resulting prices are 
therefore the Locational Marginal Prices [48 1. This implies 



that a generating unit injecting power at a given node is paid 
the locational marginal price corresponding to that node; and 
conversely, a demand receiving power from a given node pays 
the locational marginal price corresponding to that node. The 
most dominant network constraints are due to line capacity 
limits and network losses. Technical constraints cause the 
power flow through any line to be limited. The power flow is 
said to be congested when it approaches its maximum limit. 
This constraint is explicitly included in our model below. The 
second constraint is due to losses, most of which are due to the 
heat loss in the power lines. For ease of exposition, such ohmic 
losses are not modeled in this paper. The standard market- 
clearing procedure is based on Social Welfare [46 1 . Denoting 
this as Sw, it is defined as 

N D j N Gi 

^=EE ^DJkPDjk "EE X G, b P G*b (25) 
j£D q k=l ieG s b=l 

where the first and second term denote the revenue due 
to surpluses stemming from bids from GenCo and ConCo, 
respectively. The market-clearing procedure is then given by 



Maximize S 



(26) 



subject to 

N G i N Dj 

EE^-EE^w- 

iee 6=1 je-d ft=i 

E B nm [& n ~ <5m] = 0; p^n G N 



(27a) 



m GO 



B nm [S n - 5 m ] < P™ x ; 7»mVn € N;Vm e £1 (27b) 



The constraints ( |27a| i and ( |27b| i are due to power balance and 
capacity limits, respectively. It can be seen that the associated 
Lagrange multipliers, p„ and 7„ m , are indicated in each 
constraint. 

The underlying optimization problem of the ISO can be 



therefore defined as the optimization of (26i subject to con 



straints (27ai and (27bi. This problem can be restated as 
the solutions of the generation power blocks levels Pcib, the 
demand power blocks levels Pojk, the voltage angle S n , and 
dual variables p n and j nm such that 



P G ib{\ 



Gib ~ Pn{i)) 



OVi G Gf]b=l 



Pn4k(Pr,.a-\ - ^ojk) = Vj <G D q ; k = 1, . 

/ ■-' nm \_Pn Pm\ ~t~ 



Djk{PnU) ~ *D]k) 
Sn( E P" 



E 



B nm [inm ~ Imn] ) = Wl G N 



,N G i (28a) 
-,N Dj (28b) 



(28c) 



Nr. 



N 



D.1 



( - E E p « b + E E p ^k+ 

iee 6=1 jetf fc=i 

E - g m] ) = Vn G N; p n > (28d) 

Inm (P;Z X - Bnm [5 n ~ 5 m ] ) = Vn G N; Vm G n (28e) 

The decision variables of this problem are the amounts of 
power to be generated by each generating unit i in each 



block b, i.e., Pab', the amounts of power to be consumed by 
each demand j in each block k, i.e., Pnjk an d the locational 
marginal prices, p n uy It should be noted that (28di is arrived 
at using the balance equation, which implies that the term in 
the corresponding parenthesis is zero, and by multiplying this 
term by p n , which is positive. 

In summary, the market model that we consider in this 
paper include optimization goals of GenCo discussed in sec- 
tion II-A. and formulated in eqs. (lli-(12ci 
goals of ConCo, (fT5}-( |16c| i 
finally optimization of the ISO discussed in section II-D 



optimization 
and formulated in II-B, and 



and formulated in eqs (25l-(27bi. The combined problem 
statement is the determination of an equilibrium point which 
collectively optimizes all the above three market contributors. 
The corresponding equilibrium point is therefore defined as 
that which satisfies the following criteria: 

1) maximum profit for every individual generating com- 
pany 

2) maximum utility for every individual consumer 

3) maximum net social welfare for the ISO. 

IV. Market Equilibrium 
We now discuss the overall market equilibrium. In the 



GenCo and ConCo optimization problems, given in ( 1 1 i-( 12c I 



and ( 15 H 16c I, respectively, the Locational Marginal Prices 
appear as inputs. On the other hand, these LMPs appear in the 
optimization problem of the ISO in p5j)-( |27b} as dual prices 
corresponding to the balance equations, causing a tight link 
between the three families of optimization problems, thereby 
making the overall problem nontrivial. 

We note that each of these three sets of optimization 
problems are linear programming problems. Thus, the Karush- 
Kuhn-Tucker (KKT) optimality conditions are both necessary 
and sufficient for describing the optimal solutions. The op- 
timality conditions of the three sets of problems result in 
a Linear Complementarity Problem (LCP) which correspond 
to the market equilibrium for all i £ G/, b = 1, . . . , Na, 
j £ D q , k = 1, . . . , Njjj, n £ N, and meSlas follows: 



Pdb{ - Pn{i) + A 



+ OLi + (f>ib) = 



b=l 
f pmax 



Gib 



Vib {^Gib ~ P Gib 
PD]k{Pn{j) - A 

i>3 " 



Gi ) 



= 



D 



k=l 



r>max\ 
^Djk ) 



Pn(i) 

PDjk {Pn(j) - 

s n ( B nm [p 



^Djk) 



^ ^ Brim 



In 



J)=o 



(29a) 

(29b) 

(29c) 
(29d) 

(29e) 

(29f) 
(29g) 
(29h) 

(29i) 



N Gi N Dj 
i<=6> h=l j6i? fc=l 



B nm [S n - S m ] ) = 



7n 



/ pmax r> [r _ r 1\ r, 



(29j) 



(29k) 



The above problem can be compactly stated as the solution 
of the following LCP: Find a vector x* € R" that solves the 
following constraints: 



X T (Mx + q) = 
x > 

Mx+q > 



(30) 



where x 6 R nxl is the variables vector and square matrix 
M £ jj nx ™ anc i q g R nxl contains operating costs and bids of 
generators, utilities and bids of consumers, maximum and min- 
imum limit of generators and consumers and also maximum 
thermal limit of transmission lines (See Nomenclature for all 
definitions). The corresponding solution x* of the underlying 
LCP problem, which determines the market equilibrium and 
is dependent on M and q, is denoted as 



LCP(M, q). 



(31) 



Definition 10: A matrix M £ R nxn i s called a P-matrix if 
its all principal minors are positive. 

We now state and prove the theorem that discusses the LCP 
solution. 

Theorem 3: In LCP problem defined in pO) , matrix M is a 
P-matrix if and only if the LCP(M, q) has a unique solution 
for any q £ R n . Moreover, if M is a P-matrix then there is 
a neighborhood M. of M, such that all matrices in M. are 
P-matrix. 

Proof: see El. ■ 
An additional result that will prove to be useful for our 
following discussions concerns the dependence of x(M, q) on 
the market parameters. This is characterized in the following 
theorem. 

Theorem 4: |55 1 Matrix A is a P-matrix if and only if 
(I — D + DA) is nonsingular for any diagonal matrix D = 
diag(d) with < di < 1, i = 1, 2, . . . , n. Furthermore, if 
x(A, q) and x(B,p) are the solution of the corresponding 
LCP(A,q) and LCP(B,p) respectively, we have: 

|| x(B,P)-x(A,q) || 



x(B,P) 



<P(A) || (A-B)x(B,p)+q-p\\ 

(32) 



where 



(3(A) = max dem \\(I -D + DA)~ l D\\. (33) 

In the following theorem, the connection between market 
equilibrium and Pure Strategy Nash Equilibrium is presented. 

Theorem 5: Assume that all GenCos and ConCos are ra- 
tional players. If M is a P-matrix, then equilibrium of the 
market denoted as x* = LCP(M, q) exists and is identical to 
a unique Pure Strategy Nash Equilibrium. 



Proof: We prove this proposition by using contradictions. 
That is, assume that GenCos and ConCos are rational, M is 
a P-matrix and x* = LCP(M, q) is not a Pure Strategy Nash 
Equilibrium. Since M is a P-matrix, then Theorem [3] implies 
that x* = LCP(M, q) exists and is the unique maximizer of 



a constant r] as 



the GenCo problem denoted in (Hi - ( 12b i, ConCo problem 



denoted by ( [15} - ( |16b| i, and ISO problem denoted by p6 
d27b"l ). That is 



u^x^x*^) < u i (xt,x M L i ) Vxi e Si 



(34) 



Since GenCos and ConCos are rational, the inequalities in 
( fT4| > and ( fT8| i hold. From Theorem 2, it follows that the payoff 
functions Uj) are diagonally strictly concave. Theorem 

[2] in turn, implies that the game between GenCo, ConCo, and 
ISO has a unique Pure Strategy Nash equilibrium denoted as 
xn*. Using the Nash Equilibrium definition in it follows 
that 

Ui(xi,x%) < Ui(x* N .,x* N _.) Vxj € /Sj. (35) 

Since we assume that x* = LCP(M, q) is not a Pure Strategy 
Nash Equilibrium, it implies that both x* and x* N maximize the 
payoffs which contradicts the uniqueness of x* = LCP(M, q) 
and therefore completes the proof. ■ 

Remark 1: Theorem [5] provides conditions under which the 
market has a unique equilibrium, i.e. no degeneracy can occur 
lfl4l in the market if M is a P-matrix. 

We now proceed to the effect of uncertainties and curtail- 
ment factor on the market equilibrium. 

V. Perturbation Analysis of Market Equilibrium 

Using the results in the previous section, we now present 
a perturbation bound for the equilibrium of the electrcial 
market. As we discussed above, the market equilibrium is 
the solution of LCP(M ,q) which is denoted as x* and has 
been explicitly defined in Eq. p0| ). As we discussed in Section 
II. C, the wind forecast error defined in (20 1, parametrized by 
Agw, and demand response defined in ( 23 i, parametrized by 
KDjk, are considered as two sources of perturbations. These 
two components of uncertainty can affect M and q which 
explicitly in our analysis are denoted by M +AM and q+Aq, 
respectively. Therefore the underlying LCP problem under 
parametric perturbation is altered as 







x T ((M + AM)x + (q + Aq)) 
x>0 

(M + AM)x + (q + Aq)>0 
This in turn leads to a corresponding equilibrium 

x* A = LCP(M + AM, q + Aq). 



(36) 



(37) 



In the rest of the paper we are looking for the maximum shift 
of the equilibrium. The following definitions are useful. 

Definition 11: Define non-dimensional perturbation param- 
eters cm and e q as 



||AM||<e M || M||, 
l|A«||<e, || q ||, 



(38) 
(39) 



r l = e M p{M) || M||, 



and a set A4 as 



M : {AM|/3(M)||AM|| < 77} 



(40) 



(41) 



We now quantify the relation between x* A in ( 37 1 and x 



in ( 3 1 1, in the following theorem. 



Theorem 6: If the nominal market ( 30 1 has unique solution 
and rj < 1, then the perturbed market in <j36j has a unique 
solution x* A and satisfies the following inequality 

|| x* - x\ || . 2e 



x* || 1 — 1] 



(42) 



Proof: If the nominal electrical market ( |30| > has a 
unique solution according to Theorem 1, the matrix M is 
P-matrix. This in turn, according to Theorem 2, implies 
that (I — D + DM) is nonsingular for any diagonal matrix 
D = diag(di) with < di < 1, i = 1, 2, . . . ,n. We have the 
following equality 

(I - D + D(M + AM)) = (I - D + DM) (I + M AM) 

(43) 

where 

M = (I - D + DM)- X D. (44) 
Noting the definition of /3(M), it follows that 

||M AM|| < (3{M)\\AM\ \ < r\ VAM e M (45) 

Since rj < 1 according to theorem 3, it follows that I+M AM 
is nonsingular for all AM g A4. We therefore conclude from 
(43} that (I-D + D{M + AM)) is a P-matrix and therefore 
the market in (36 1 has a unique solution for all AM G Ai 
and r) < 1. 



We now show that x* A satisfies (42 1. Using (32 1 in Theorem 
2, we have that: 



II x * t < P(M + AM)\\AMx* + Aq\\ 
We rewrite the argument of /3(M + AM) as 
(I - D + D(M + AM)y X D = 
(I + M (AM)) _1 (J - D + DM)^D 
Using Appendix 1 in 0, we have: 

,-!„ . 1 . 1 



(46) 



||(7 + M (AM))" 



< 



< 



1-/?(M)||AM|| - 1-r/ 
Taking norms on both sides of ( |4"7} , we obtain that: 



P{M + AM) < 



1 



P(M) 



Using (46 1 implies that: 

II x* 



< 



1 



- p(M) 1 1 AM + Ag| | 



(47) 



(48) 



(49) 



(50) 



II x " II 1 ~ V 

Considering the definition of AM and Aq, we obtain, in turn, 

|| x* - x\ || . 2e 



< 



1 



(51) 



where e = max{eM || M \\,e q \\ q ||} which is the desired 
bound. ■ 

Remark 2: Defining /ifAga, KDjk) f°r a U I <= G w , b = 
l,...,N Gi ,jeD q ,k=l,. 



TABLE I 

Cost functions data of Generators 



fJ,(A G i b , KDjk) 



,N Dj as 
2r 



l-rj 



(52) 



Theorem [6] implies that the uncertainty in market can lead 
to a maximum shift by the equilibrium of an amount 
fi(Acib, KDjk)- As this function is nonlinear, determining 
an analytical relationship between (i(AGlb> K Djk)> Aaib an d 
KDjk is exceedingly difficult. As will be shown in the next 
section, simulation studies show that as KDjk increases, ji 
decreases. This in turn brings the perturbed equilibrium closer 
to the nominal equilibrium. 

A. Game theoretic Interpretations 

We now evaluate the perturbed market using tools from 
game theory. Consider two strategic games defined by two 
profiles of utility functions as G = (X, (S'i)igi, (tt,),gi) and 
G = (I, (Si)i e x, (uj)jgi). If s* is a Nash equilibrium of G, 
then s* need not be a Pure Strategy Nash Equilibrium of G. 
The equilibria of the strategic games G and G may be far apart, 
even if (uj)igx and (ttj)i£z ai" e ver Y close to each other. 

Definition 12: Given e > 0, a Pure Strategy Nash Equilib- 
rium s* E S is called an e — equilibrium if for all i e I and 



Ui(si,s*_ J < Ui(s*,s*J + eVsi £ Si 



(53) 



Obviously in Definition [12] when e = 0, the e — equilibrium 
is a Pure Strategy Nash Equilibrium in the sense of Definition 

m 

The following definition formally defines the closeness of 
two strategic form games. 

Definition 13: Let G = (X, (Si)i^x, (uiji^x) and G — 
(I, (Si)iez> (wi)iez) be two strategic form games, then G is 
a a — approximation of G if for all i el, 



|u,j(s) — Ui(s)\ < a Vs € 5. 



(54) 



The next proposition relates the e — equilibrium of close games 
as defined in [13] 

Theorem 7: If G = (1, (Si)igz, (u»)iez) is an a - 
approximation to G = (I, (<5j)j 6 z ; (ui)iex) an d s* € S is 
an equilibrium of G, then s* e 5 is a (2a) — equilibrium of 
G. 

Proof: For alH € I and all € S*i, we can write 
Ui(Sj, S*J - Ui(s*,S*j) = Ui(Si, s*J - ^(s^s*,)-!- 

(55) 

Since s* is a Pure Nash Equilibrium of G, then Wj(sj, sl j — 
Ui(s*,s*J < 0. This in turn follows that 



Name 


Block 1 

Sizef^™') 


PriceWMW) 


Block 2 

Size< Mlv ) 


p rice ($/MMO 




12.5 


240 


40 


900 




12.7 


300 


40 


950 




12.2 


240 


80 1200 




12.2 


240 


80 1200 


Pgr.;, 


10 


40 


20 


80 



Using Definitions [13] and [12] it follows that 

Ui(si,s*_ t ) - Ui(s* , s*_i) < 2a 



(57) 

equilibrium of 



which in turn implies that s* g S is an (2a) 
G. ■ 

Proposition 1: Assume that the perturbed market denoted 
by the strategic game G under uncertainties Agn and KDjk 
is an a — approximation of the nominal market denoted by G. 
If x* A is an equilibrium of the perturbed market G, then x* A 
is an (2a) — equilibrium of the nominal market G. 

Remark 3: If the wind forecast error Acib increases, from 
Definition [13] it follows that a will increase as the relative 
payoffs of the corresponding generators will have an increased 
error. Proposition[T]implies that in such a case, the equilibrium 
of the perturbed game G is correspondingly far away from 
the nominal market G. It should be noted that as the forecast 
error increases, the corresponding cost of deploying ancillary 
services increases as well. If a Demand Response program is in 
place, this cost increase is conveyed to the consumer, leading 
to a decrease in the load quantified by the demand curtail factor 
KDjk as in (23 1. Similar to our observation in Remark [2] it 
should be noted that the exact impact of an increasing KDjk 
on a is difficult to quantify due to its nonlinearity. However, 
simulation studies show, as discussed in Section |VI| that as 
KDjk increases, a decreases. This in turn brings the perturbed 
equilibrium closer to the nominal equilibrium. 

VI. Case Study 

An IEEE 30-bus case is used for simulation studies, whose 
interconnections are shown in Figure [T] The size and price 
of each block of each GenCo are shown in Table [I] and for 
the sake of simplicity we assume that each GenCo bids in 
its marginal cost. The size of minimum power requirement of 
each demand and the size and price of each block of each 
demand are shown in Table [II] We assume that only part of 
the overall loads are participating in the market, and the fixed 
loads in other buses are shown in Table [III] The reactance B nm 
of the line connecting bus n and bus m can be found in Table 
IV The transmission capacity limit of all lines is chosen to 
be 100MW. The line parameters are perunit with three-phase 
base of 230 kV and 10 MVA. 

A. Nominal Market Equilibrium 

Table [V] provides equilibrium results concerning generator 
output, revenue, and profit for the electric power market if 
no wind shortfalls are imposed to the system. These results 



(56) are obtained by directly solving LCP problem (29ai-(29ki 



TABLE II 

Cost functions data of Consumers 



Name 


Block 1 


Price( $ /MW0 


Block 2 

Size< MW/ ) 


Price ($/A/wo 


Pd 7 


11 


150 


19 


600 




10 


100 


21 


600 


Pdn 


12 


170 


22 


600 


Pdo 


12 


150 


24 


600 


Pd-26 


15.5 


150 


21 


600 


Pd 2 . 


9.5 


150 


22 


600 



TABLE III 
Fixed load data 



Name 


Demand 


Name 


Demand 


Pd, 


22.31 


Pd^ 


7.20 


P di 


8.83 


Pd le 


4.06 


Pd» 


13.95 


Pd 17 


10.46 


p d 10 


6.74 


Pd 18 


3.72 


P dl2 


13.01 


Pd, a 


11.04 


Pd-20 


2.55 


Pd 2l 


3.39 


Pd23 


22.31 


Pd 2i 


10.11 


Pd 2 e 


4.06 


Pd 2 g 


2.78 




Table [VI] shows the power consumed and the corresponding 
demand payments. Since we ignore losses in this study, the 
difference between the total generator power output and the 
total power consumed, and are equal to zero and both are 
221 whereas 15l MW of the load is a fixed demand as 



can be seen in Table III The Locational Marginal Price in all 
of the buses is 23.62F^ 4Wh due to the fact that transmission 
lines are not congested. 



Fig. 1. IEEE 30-bus case study 

TABLE V 

Results of market equilibrium for GenCo, no wind uncertainty 



Name 


Power output 1 ™' 


Revenues*/ h 


Cost $ / h 


Total Profit*/' 1 




30 


708.6 


660 


48.6 


p*2 


11.64 


274.93 


274.76 


0.18 




80 


1889.6 


1000 


889.6 


Pgn 


80 


1889.6 


1000 


889.6 


Pg 13 


20 


472.4 


8 


464.4 



TABLE IV 

Transmission Lines Reactance; Line reactances are perunited 
based on lOOMiy 



Connected 




Connected 




Bus 


Reactance 


Bus 


Reactance 


Bus 


Bus 


(x b _k) 


Bus 


Bus 


(Xb_k) 


b 


k 




b 


k 




1 


2 


0.0575 


1 


3 


0.1652 


2 


4 


0.1737 


3 


4 


0.0379 


2 


5 


0.1983 


2 


6 


0.1763 


4 


6 


0.0414 


5 


7 


0.116 


6 


7 


0.082 


6 


8 


0.042 


6 


9 


0.208 


6 


10 


0.556 


9 


11 


0.208 


9 


10 


0.11 


4 


12 


0.256 


12 


13 


0.14 


12 


14 


0.2559 


12 


15 


0.1304 


12 


16 


0.1987 


14 


15 


0.1997 


16 


17 


0.1923 


15 


18 


0.2185 


18 


19 


0.1292 


19 


20 


0.068 


10 


20 


0.209 


10 


17 


0.0845 


10 


21 


0.0749 


10 


22 


0.1499 


21 


22 


0.0236 


15 


23 


0.202 


22 


24 


0.179 


23 


24 


0.27 


24 


25 


0.3292 


25 


26 


0.38 


25 


27 


0.2087 


28 


27 


0.396 


27 


29 


0.4153 


27 


30 


0.6027 


29 


30 


0.4533 


8 


28 


0.2 


6 


28 


0.0599 









B. Perturbed Market Equilibrium with Wind Uncertainty 

We assumed that GenCo at bus 13 in figure [2] is wind 
based and is subjected to an uncertainty Agi3- For the sake 
of simplicity in simulations we assume that uncertainty for all 
blocks are the same and characterized by A<3 13 . The resulting 
perturbed market, discussed in Section |VJ was simulated. 
Figure [2] shows the corresponding market equilibrium shift 
[i(Acib: i^Djk) as define in d52l for different wind penetration. 

As can be seen in Fig |2J with increasing Agi3, higher 
equilibrium shift [i results. As wind penetration x w % defined 



in (22 1 increases, the corresponding market equilibrium shift 



TABLE VI 

Results of market equilibrium for ConCo, no wind uncertainty 



Name 


Power consumed 


Total Payment % / h 


P*r 


11 


259.82 


Pdu, 


12 


283.44 


Pd:m 


10 


236.2 


Pd, 


15.5 


366.1 


Pd26 


9.5 


224.4 


P d27 


12.0 


283.44 




Fig. 2. Market 
penetration 



equilibrium shift as a function of uncertainty and wind 




Fig. 3. Market equilibrium shift as a function of uncertainty and demand 
curtailment factor 



deviation is higher. 

C. Measure of Uncertainty incorporating Wind Uncertainty 
With Demand Response 

We now introduce a perturbation Kdis into the picture 
to represent demand response which correspond to introduce 
Real Time Pricing (RTP) at the consumer at bus 15. We now 
numerically evaluate the effect of this perturbation parameter. 
In particular, we evaluate the equilibrium shift fi(Aaib, KDjk)> 
for a range of npi^ at different level of wind forecast error in 
Agi3. Figure [5] shows the corresponding results, which clearly 
demonstrate that RTP can reduce the equilibrium shift and 
therefore the cost of uncertainty in the presence of a wind 
forecast error. 

VII. Concluding Remarks 

The current energy crisis has created an urgent need in 
integrating renewable energy resources into the power grid. 



The latter in turn can introduce intermittency and uncertainty 
into the picture, thereby introducing a prohibitive integration 
cost. In this paper, we introduce an analytical framework 
to evaluate this cost. We begin with an overall model of 
the energy market including Generators Company, Consumers 
Company as well as Independent System Operator (ISO). 
We analyze the underlying market equilibrium using game 
theory and established sufficient conditions are derived for the 
existence of a unique Pure Nash Equilibrium for the nominal 
market, then the effect of uncertainty due to renewable energy 
on this equilibrium using perturbation analysis. The perturbed 
market is analyzed using the concept of closeness of two 
strategic games and the equilibria of close games using the 
notion of a — approximation and e — equilibrium. This analysis 
is used to quantify the effect of uncertainty of RERs and its 
possible mitigation using Demand Response in the form of real 
time pricing, we quantify its effect on the market equilibrium 
using a parameter denoted as Curtailment Factor. We show 
that the equilibrium shift caused by the RER uncertainty is 
mitigated in the presence of RTP due to the Curtailment Factor. 
Numerical results are included that validate the theoretical 
results, using an IEEE 30-bus network. 
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